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A universal algorithm for a deterministic preparation of arbitrary three-mode 
bosonic states is introduced. In particular, we consider preparation of entangled 
quantum states of a vibrational motion of an ion confined in a 3D trapping poten- 
tial. The target states are established after a proper sequence of laser stimulated 
Raman transitions. Stability of the algorithm with respect to a technical noise 
is discussed and the distance (fidelity) of outputs with respect to target states is 
studied. 



1. Introduction 

The deterministic control of preparation of genuine quantum states is one of the 
desired goals of quantum physics. Recent advances in atomic physics (trapped ions [1, 
2]) and quantum optics (cavity QED [3]) have enabled an effective control of dynamics 
of microscopic quantum systems. These systems can serve for tests of fundamental 
concepts of quantum mechanics and quantum measurement theory. Other applications 
can be related to quantum computing [4, 5] and information processing [6]. 

A great effort has been devoted to the task of preparation of various nonclassi- 
cal states (e.g., Schrodinger cats, Fock and squeezed states) [7]. Formal introduction of 
some nonclassical states (e.g., squeezed states) has been based on the Perelomov scheme, 
i.e., under the action of an appropriate unitary operator on a reference state. This most 
simple way of quantum-state preparation can be represented as a time evolution of an 
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input state governed by a proper Hamiltonian. However, physical implementation of a 
required Hamiltonian can face insurmountable difficulties. Quantum-state-engineering 
methods introduced firstly in the cavity QED (micromasers) have overcome these re- 
strictions by performing conditional measurements on atoms which result in a projection 
of the cavity field on desired target states [8, 9]. However, due to a specific selection of 
quantum trajectories these methods suffer from low probabilities of desired outputs. 

The process of deterministic preparation of an arbitrary one-mode state of electro- 
magnetic field has been described recently by Law and Eberly [10]. In their approach 
a proper sequence of unitary transformations associated with a sequential switching of 
two interaction channels transforms the vacuum state of the field to prescribed tar- 
get wave-function. The experimental preparation of quantum states of 1-D vibrational 
motion of trapped ions has utilized switching between interaction channels via tuning 
lasers either to electronic transitions or appropriate vibrational sidebands [1]. Possible 
generalizations of the 1-D quantum state synthesis to 2-D have been discussed recently 
[11, 12, 13]. 

In this paper we consider a deterministic preparation of multi-mode bosonic states. 
In particular, we consider synthesis of quantum states of a vibrational motion of an ion 
confined in a 3D trapping potential. The main motivation for utilizing trapped ions is 
a possibility to neglect decoherence effects due to dissipations. An effective control of 
coupling between external and internal degrees of freedom via a proper laser irradiation 
is also advantageous. 

2. Quantum state synthesis 

Let us consider a preparation of an arbitrary three-mode bosonic state 

M x M y M z 

\$ 'target) = ^2 ^2 ^Qn x ,n y ,n z \n x ,n y ,n z ). (1) 
n x —0 n y —0 n z — 

In particular, the three modes correspond to vibrational modes of a quantized center- 
of-mass motion of an ion confined in the 3D trapping potential. The key idea of 
the quantum state synthesis is a splitting of the whole Hilbcrt space H V ib of three- 
mode vibrational states into a hierarchy of two-dimensional subspaces. This enables 
us to transfer populations between pairs of component Fock states (see below). How- 
ever, this dynamical "ordering" requires a coupling of vibrational degrees of freedom 
to an adjoint system. For the trapped ion it is natural to assume the coupling be- 
tween external (vibrational) and internal (electronic) degrees of freedom. For purposes 
of our quantum-state-synthesis algorithm we utilize four internal electronic levels \i) 
(i = a,b,c,d) spanning the Hilbert space Tii n - In general, during the preparation the 
vibrational and internal states become entangled and the whole state vector can be 
written as a superposition of "component" states 

M x M y M z 

!*(*)> = EEE E Qn x ,n y ,n z ;i(t)K,n y ,n z )^\i) (2) 
n x =0 n y =0 n z =0 i=a,b,c,d 
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Fig. 1. The scheme of the de-evolution procedure (3): a recursive transfer of excitations to a 
subspace with a number of trap quanta decreased by one. 



Jmax J J—n x 

E E E 

J—0 n x —0 n v —0 i- 



a.b.c.d 



-n y \i {t) \ n x 



Tly , J 



<8> |i>. 



The first summation in the second expression is taken over subspaces 7ij of the vi- 
brational Hilbert space TL V ib- The subspaces TLj are spanned with component states 



\Tl x , Tl 1 



, n z ) with a fixed total number of trap quanta J = n x + n y + n z 



Firstly we solve for the given target state (1) an inverse task, i.e. we "de-evolve" 
it towards the reference vacuum state via a proper sequence of "elementary" unitary 
operations which are associated with particular laser stimulated Raman processes (see 
below). The "de-evolution" into vacuum can be schematically written as (see Fig. 1) 



|0, 0, 0> <8> |a) 



fjabc 
U 



fjbed 
J-1,J U J 



l U J 



be 



J=l 



I* 



target 



) ® \a). 



(3) 



Here superscripts and subscripts indicate electronic levels and vibrational states from 
subspace Tij which are involved in a controlled "de-evolution" of the state vector (2). 
Namely, Uj bc transfers populations of components states \n x ,n y ,n y ) <g> \i) with n x + 
n y + n y = J an( i i = a, b,c to a, single component state \J, 0,0) ® |a) which then 
enters the state vector (2) on the subspace under consideration. Next U f_ x completely 
transfers populations of components states \n x , n y , n y )®\i) with n x +n y +n y = J— 1 and 
i = b, c, d to the state | J— 1, 0, 0)® \b). After that Uj b _ 1 j transfers the actual population 
of the components state |J, 0,0) ® \a) to |J — 1,0, 0) (g> \b). As result, the sequence 
Uj b _ 1 jUj cd : 1 Uj bc outlined in Fig. 1 cancels the contribution of component states from 
the Hilbert subspace Hj (E> Hi n into the state vector (2). Recursive application of this 
sequence on subspaces with decreasing J evolves the initial target state into the vacuum. 

Note that the solution of the inverse task (3) offers immediately also the prescription 
for preparation of the target state from the vacuum as 



target/ 



\a} = 



?7 tabc |0,0,0) 



(4) 
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In what follows we describe in details how to "construct" unitary transformations ap- 
pearing in Eq.(3) from a set of "elementary" transformations (time evolution opera- 
tors) associated with laser stimulated Raman processes for a trapped ion. Therefore 
the preparation sequence (4) corresponds to the action of "elementary" transformations 
applied in the opposite order and with phases of lasers shifted globally by ±tt. 

3. Manipulation of trapped ions: laser-stimulated Raman processes 

In this section we briefly introduce "tools" necessary for our manipulations with 
vibrational states of trapped ions. For quantum state synthesis of 3D motional states 
of a trapped ion we utilize laser stimulated Raman processes [14]. Let us consider a 
trapped ion confined in a 3D harmonic potential characterized by trap frequencies v q in 
orthogonal directions q = x,y, z. The ion is irradiated along two axes (y and z for con- 
creteness) by two external laser fields with frcquencie and wave vectors k y , hz ; 

respectively. The laser fields stimulate Raman transitions between two internal energy 
levels | a) and \b) via an auxiliary electronic level which is far off the resonance. After 
the standard dipole and rotating-wave approximations (RWA at laser frequencies) the 
adiabatic elimination of the auxiliary off-resonant level leads to the effective interaction 
Hamiltonian [14] 

H{ eff) = g*e- i ^-^ t D y (-ie y )D z (ie z )\b)(a\+h.c. (5) 

Here D q (ie q ) = e ie i^ a l +a ^ = Q lk <iQ [ s the displacement operator; a g ,aj are the creation 
and annihilation operators of the vibrational mode q. The corresponding Lamb-Dicke 
parameter e q is defined as e 2 q = % k 2 /{2mhv q ). The effective interaction constant g ~ 
£*£ Z /A is proportional to complex laser amplitudes £ q and inversely proportional to 
detuning A of lasers from the auxiliary level. The resonant terms in the Taylor series 
expansion of Eq. (5) which contribute dominantly to the resulting effective Hamiltonian 
can be selected by an appropriate choice of laser frequencies. If off-resonant processes 
are oscillating with sufficiently high frequencies they can be eliminated applying the 
second RWA at trap frequencies. In particular, tuning lasers to the first vibrational 
sidebands in order to adjust resonance conditions ui\ = ui z — u) y = LUb — u a + v y — v z we 
find that the retained resonant terms lead to the interaction Hamiltonian [14] 

H x = glalT ey [n y }AAnz}a z \b)(a\e- iu;it + h.c., (6) 

where 

T tq [a\a q ] = exp(-- e 2)£^_J^at^. (7) 

The generalized Rabi frequency £ln],n z of the transition between the states \n x , n y , n z )® 

\a) and \n x ,n y + l,n z -l)®\b) is given as e - 2 ( e y)^-i( e z)/v(% + l) n * where 
L\ is the associated Laguerre polynomial. In the Lamb-Dicke limit of e q <C 1 it reads 
Clny,n z — * \/ (% + l)n z as T Eq is close to the unity operator. 
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With the proper laser tunings we can design also other interaction Hamiltonians 
required for the quantum state synthesis. Namely we consider counterpropagating laser 
beams in the x direction which classically drive the electronic transition |o) <-> \b) via a 
resonant laser-stimulated Raman process. The resonance condition for laser frequencies 
uu 2 = u x2 — bj x \ = uu}> — w a leads to the interaction Hamiltonian 

H 2 = g^\b){a\e-^ t +g 2 P < \a){b\e^ t . (8) 

The effective Lamb-Dicke parameter e' x with k' = k x2 — k x i is much smaller compared 
with e x and thus for simplicity it is assumed to be close to zero. Consequently can be 
approximated by the unity operator. The analogous process with the laser frequencies 
such that u)g = lo x2 — u> x i = uJb — uj a — u x is described by the one-mode interaction 
Hamiltonian H3 

H 9 = fl ^ e; a I |6>(o|e- <w » t + flfl^ot loX&le*-* (9) 
which is known as the nonlinear Jaynes-Cummings model [15]. 



4. "De-evolution" procedure as sequence of "elementary" transformations 

Let us now assume that at some stage of the "de-evolution" procedure we operate 
on the subspace Tij ® Habc and contributions of component states from subspaces with 
higher total number of trap quanta J as well as from the subspace Hj <8> Hd into the 
state vector (2) have been already canceled. (The subspaces of Hi n are spanned by 
indicated electronic levels.) Firstly, we explicitly construct the operator Uj bc from 
Eq.(4) which transfers completely populations of component states from the subspace 
under consideration to the single component state | J, 0, 0) £g> \a). The operator Uj bc can 
be expressed as 



.7-1 



jjabc 
Aj,n x 



(5) 

\n x , J— n x ,0;b) ' 



n x =0 



u {2) u {1) u (2) u {1) 

|n x , J— n x ,0;b) \n x ,J— 1— n x ,l;a) '** \n x ,l, J— 1 — n x ;b) 1^,0, J— n x ;a} ' 



(4) 



(3) 



\n x ,J—n x ,0;c) \n x , J— 1 — n x ,l;b) 



(4) 



(3) 



\n x ,l y J— 1— n x ;c) |n x ,0, J— n x ;b) 

tip) 



(10) 



(4) 

\n x ,0, J — n x ;c) ' 



Here an "elementary" unitary transformation n n = exp(—iTH p ) is equal to the 

evolution operator corresponding to a particular interaction Hamiltonian H p applied for 
time r. Subscript indicates a component state \n x , n y , n z )®\i) the population of which is 
completely transferred to another component state according to the applied interaction 
"channel" p. 

Actions of the "elementary" transformations which form the operator Uj bc are out- 
lined in Fig. 2a. In particular, under the action of J _ n O . o ^ [see Eq.(10)] the 
complete transfer of the population of the component state \n x , 0, J — n x ) ® \a) to the 
component state \n x , 1, J — 1 — n x ) <S> \b) is achieved providing that the interaction time 
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Fig. 2. The scheme of the de-evolution procedure (3). (a) Transfer of populations within the 
subspace with a given number of trap quanta J under action of U] bc . (b) The "neighbor- 
ing" subspace with the number of trap quanta J — 1 under action of Uy^. Labels of arrows 
indicate the "direction" of action of the elementary transformations C/' p '. 



t and parameters |<7i|, 0\ of the Hamiltonian H\ (6) are chosen to fulfill for n y — the 
condition (n z = J — n x — n y ): 

ie- iei Q„ x ,„ B ,„ i; aCos[| 5l |Ttt^ (11) 

Next, the transfer of the population from the component state \n x , 1, J — 1 — n x ) <g> \b) 
to the component state \n x , 1, J — 1 — n x ) ® |a) is achieved under the action of the 
interaction Hamiltonian H2 (8). The following choice of the interaction parameters 
\g2\T and 82 is required for n y = 1 

^ ie2 Qn x ,n y ,J-n x -n y ;bCOs(\g2\r) + Qn x ,n y ,J-n x -n y ;a sin(|0 2 |T) = 0. (12) 

The repeated sequence of these "elementary" transformations forms the operator Aj, n:c 
[see Eq.(10)] which transfers the population of the states \n x ,n y , J — % — n x ) (g> \i) with 
n y = 0, . . . , J — n x and i = a, b to the component state \n x , J — n x , 0) (g) |a). 

Further, under the action of the operator Bj jTlx+ i populations of states \n x + 
1, n y , J — n y — (n x + 1)) ® \i) with n y = 0, . . . , J — (n x + 1) and z = b,c are trans- 
ferred to the component state \n x + l,J— (n x + 1), 0) ® It is seen from Eq.(10) that 
the form of -B/.^+i (acting on "neighboring lines" in Fig. 2a with number of x-quanta 
increased by one) is quite analogous to ^4j,n x but instead of levels \a), \ b) the electronic 
levels 1 6), |c) are involved in the transfer of populations. The corresponding Hamiltoni- 
ans i?3, H4 and the transfer conditions are obtained from Eqs.(6),(ll) and (8), (12) by 
the substitutions: \a) — » \b) and \b) — > |c). 



Quantum synthesis of 3D vibrational states of trapped ions 



277 



As the next step of the "de-evolution" sequence (10) the operator Cj t n x — U^ n j_ n Q . b ^ 
is applied to the state vector (2). It transfers the whole population of \n x , J—n x ,0)($i\a} 
to |na:+l, J— (n x +l), 0)®|6). The corresponding interaction Hamiltonian is obtained 
from Hi after replacements y — > x, z — > y. In other words, two lasers stimulating the 
Raman transition operate in x,y directions instead of y, z. Note that the action of 
the operator Bj tTlx+ i was necessary to avoid a transfer of population backwards to 
component states with smaller n x on the given subspace Hj ® T~L a bc- 

Altogether, the role of the operator Uj bc (10) is to transfer populations of all com- 
ponent states from Hj <8> H a bc to the only component state |J, 0,0) ® \a). The next 
important step in the "de-evolution" procedure (3) is related to the operator \J h f\ 
which transfers populations of all component states from the "neighboring" subspace 
® H-bcd to its only component state | J — 1,0, 0) (g> \b). The operator Uj cd can be 
built up from a set of "elementary" transformations in a close analogy with the unitary 
transformation Uj bc by utilizing the electronic levels i — b 7 c,d instead of i = a,b,c. 
The action of U bc \ represented by a sequence of "elementary" transformations is out- 
lined in Fig. 2b. The interaction Hamiltonians H§ and Hj associated with the indicated 
"elementary" transformations are obtained from H4 and H 3 , respectively, by the substi- 
tutions |c) — ► \d) and \b) — > |c). The Hamiltonian H 8 arises from H 5 after substitutions 
\b) -» \c) and jo) -» |6). 

The next step which connects the subspaces shown in Figs. 2a, b has been already 
sketched in Fig. 1. The operator Uj b _ 1J = C^j q - a ) associated with the interaction 
Hamiltonian (9) transfers the population of the component state |J, 0,0} <£> \a) to |J — 
1, 0, 0) (g> \b). In other words, the operator Uj b _ 1 j moves us from the subspace Hj®TL a bc 
to the subspace Hj-i ® H a bc- Note that contribution from Hj-i ® into the state 
vector (2) has been canceled by the previous action of \J h j c \- It means that we can 
apply the whole "de-evolution" procedure (3) recursively decreasing the total number 
of trap quanta. 

5. Remarks on imperfect preparation 

The deterministic preparation procedure as described above is based on the unitary 
time evolution (4) of the reference state to the target state. The solution for the 
inverse problem (3) can be found only in the absence of dissipations. However, one 
important source of noise (even for almost dissipation free dynamics of trapped ions) 
is an imperfect control of interaction constants (which depend on laser amplitudes 
and phases) and switching times required for transfer of populations of component 
states. To estimate the influence of this "technical noise" now we introduce random 
fluctuations around the ideal interaction constants and switching times defined by (4). 
In particular, we assume that fluctuations are equally distributed around the ideal values 
of interaction lengths |g p |r [e.g., solutions of Eqs.(ll,12)] within a given interval 5. Also 
phases fluctuate around their ideal values 9 P within a fixed interval 6$ . The preparation 
sequence (4) with non-ideal interaction constants and switching times then results in 
some output state \$>g) of the form (2). Its internal and external degrees of freedom 
are in general entangled. The deviation from the desired target state Y^target) ® \o) can 



278 



B.Hladky et al. 



1.0 



0.9 



□ '- 



u 

us 



0.8 



0.7 



(a) 



0.0 



A - 




□ - 


IV'>GHZ 


o 





0.01 0.02 
fluctuations 



0.03 



u 0.8 




0.6 



A - 


— l"</»diag 


□ ■ 


"~" l*>GHZ 


□ - 


P,ty»GHZ 


o ■ 


|-0> corr 


o - 


P a l*>corr 



0.0 



0.01 0.02 
fluctuations 



0.03 



Fig. 3. Fidelities of the outputs to the indicated target states in the presence of a technical noise. 

(a) Varying phase fluctuations So with fixed fluctuations S — 0.01 of the interaction lengths. 

(b) Varying fluctuations 8 of interaction lengths with fixed phase fluctuations So = 0.01. 



be quantified using an average fidelity of particular imperfect realizations \^s) to the 
target state: 



f S = ({\(*s\*tar g et)\a)\ 2 h, 



(13) 



where ((. . .)}$ stands for the average over imperfect realizations with a given range of 
fluctuations <5. We have performed simulations of noisy quantum state synthesis for 
three specific target states. We have considered the three-mode "Greenberger-Horne- 
Zcilinger" state \^ghz) — JV(\a,a,a) + \ — a,— a,— a)) with coherent components; 
three-mode correlated state \^ CO rr) = c~l"l I 2 J2 n "^fl 71 ' 71 ■> n )' w hh a = 1; and the 
superposition state \^diag) = ^ lCn=o l n ' n > n )- The cu t°ff m Fock basis is chosen 
as Jmax = 10. The corresponding average fidelities (over 100 noisy preparations) are 
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shown in Fig. 3a for an increasing interval of phase fluctuations 5$ around the ideal 
phases 9 p . The interaction lengths fluctuate around ideal values \g p \r within a fixed 
interval 6 = 0.01. The Lamb-Dicke parameters are chosen as e x = 0.3, e y = 0.1, 
e z = 0.2, e' x = 0.1. In Fig. 3b we show averaged fidelities for increasing fluctuations S 
of the interaction lengths, now with fixed phase fluctuations 6$ = 0.01. 

The figures clearly demonstrate that the preparation procedure works well even in 
the presence of relatively small fluctuations. In general, the overall error is accumulated 
by each "elementary" transformation. The total number of "elementary" transforma- 
tions scales polynomially with respect to the dimensionality of the Hilbert space. In 
our particular case the total number of operations is proportional to J^ ax - We have 
simulated the noisy preparation procedure also for higher values of the Lamb-Dicke 
parameters (up to e q = 0.4). The dependence of fidelities on amplitudes of fluctuations 
is very similar to that one shown in Fig. 3. In particular, it turns out that the fidelities 
are less sensitive to (absolute) fluctuations 5 of the interaction lengths when increasing 
the Lamb-Dicke parameters. This is due to an increase of the ideal values \g p \r of the 
interactions lengths required for the quantum-state synthesis. 

One possibility to improve fidelities of outputs is to utilize conditional measurements 
of the internal electronic levels. The preparation procedure starts and ends with the 
ion in the same internal state \a). So after preparation sequence we can verify whether 
the ion occupies the desired internal state. It can be done by driving dipole transitions 
from other internal levels to proper auxiliary levels and observing fluorescence signals. 
No signal means that the undisturbed ion is in the right electronic state \a). Wrong 
preparation sequences (with the test signal) are thrown away. In other words, outputs 
are projected on P a = \a)(a\. The improvement of fidelities due to the conditional 
selection of outputs is shown in Fig. 3b. Corresponding decrease of probabilities to 
find on the check the right final internal state is not so dramatic and state-dependent 
as known from conditional methods for micromasers. For example, the fidelity is lifted 
from 75% to 92% (64% -» 87%) for the state |* CO rr) H^ghz)) with fluctuations 5 = 0.03 
with the efficiency 82% (74%). 

6. Summary 

We have investigated preparation of entangled quantum states of the vibrational 
motion of the ion confined in the 3D trapping potential. We introduced the universal 
algorithm for deterministic preparation of quantum states which is based on a sequence 
of laser stimulated Raman transitions. In the presence of a technical noise fidelities 
of outputs with respect to desired target states can be significantly improved when 
conditional measurements on the electronic levels of the ion are performed. The method 
can be adopted also for other three-mode bosonic systems. 

Acknowledgement We thank Jason Twamley for helpful discussions. 
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